Hot Giant Loop Holography 
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We argue that there is a phase transition in the expectation value of the Polyakov loop operator 
in the large N limit of the high temperature deconfined phase of A/" = 4 Yang-Mills theory on a 
spatial . It occurs for large completely symmetric representation of the SU{N) symmetry group. 
We speculate that this transition is reflected in the D-branes which are the string theory duals of 
giant loops. 
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The Hawking-Page phase transition [ij is the col- 
lapse of hot anti-de Sitter space to an anti-de Sitter- 
Schwarzschild black hole. A beautiful picture of the holo- 
graphic dual of this transition has emerged in the con- 
text of AdS/CFT duality g-Q. On the type IIB su- 
pergravity side is black hole formation with asymptotic 
AdS5 X geometry. On the gauge theory side, we have 
the deconfinement phase transition of large N J\f ~ 4 
supersymmetric Yang-Mills theory defined on the spa- 
tial three-sphere S^. The confined phase is dual to hot 
AdS^ X whereas the deconfined phase is dual to the 
black hole geometry. 

Confining behavior of A/" = 4 Yang-Mills theory is gov- 
erned by the representation of its center symmetry and 
characterized by the expectation value of the Polyakov 
loop, which is the trace of gauge theory holonomy in pe- 
riodic Euclidean time Q Q , 



Tr U{x) = Tr Ve 



f^fdrAoir.x) 



(1) 



The loop transforms as Tr U{x) cTr U{x), where c ^ 
-g ^Yie generator of the Zn center of the SU{N) 
gauge group. Its expectation value vanishes in the low 
temperature confining phase where center symmetry is 
good, and it can be non-zero in the high temperature 
deconfined phase where the center symmetry is broken. 

The study of the Polyakov loop using effective field 
theory has a long history Q; however, it was applied to 
A/" = 4 Yang-Mills theory on a spatial only relatively 
recently Q In the weak coupling limit, all of the fields 
in A/" = 4 Yang-Mills theory have a mass gap arising from 
their conformal coupling to the curvature of the . Inte- 
grating them out to find an effective field theory in which 
one can compute the expectation value of the Polyakov 
loop is a well-defined procedure. The effective theory is 
a unitary one-matrix model with effective action S'off[C/]. 
Gauge symmetry impHes 5'off[C/] = ^off [I^C/H^'^], where 
is a unitary matrix, while center symmetry implies 
•Soffits] = •S'cff [cC^] • Further, the action is of order N'^ , 
S^s{U = 1] gjnce TV = 4 Yang-Mills theory has 

conformal symmetry, the effective action depends on the 
temperature T and the radius R through the product 
TR. It also depends on the 't Hooft coupling A = Oym-^ 
and on N which we assume is taken to infinity holding A 
fixed. At weak coupling, where A — > 0, the phase transi- 
tion is found by tuning TR to a critical value. When the 
coupling A is turned on, this phase transition is thought 



to persist and at large A to coincide with the Hawking- 
Page transition of the gravity dual. The effective field 
theory description should be reliable when TR < 1 ; how- 
ever, it is thought to have a broader applicability. We 
will assume that it can be used to discuss the deconfined 
phase, at least in the vicinity of the phase transition that 
occurs when TR ~ 1 . 

The unitary matrix model can be used to calculate the 
expectation value of the Polyakov loop operator in any 
irreducible representation R of the SU {N) gauge group. 

We shall show in the following that ^ can have interest- 
ing behavior which depends on the size and nature of the 
representation. We will consider completely symmetric 
representations Sk whose Young tableau is a single row 
with k boxes and completely antisymmetric representa- 
tions Ak whose Young tableau is a single column with 
k boxes. We shall consider large values of k so that 
remains finite as A^ — > oo. Both representations Sk and 
Ak have center charge k mod N so that ([2]) vanishes in 
the confined phase when this charge is non-zero. The 
characters can be non-zero in the deconfined phase. 

In the duality between gauge fields and strings, the 
expectation value of the Wilson loop is dual to an open 
fundamental string amplitude. This has been made pre- 
cise for the Maldacena- Wilson loop which differs from 
the Polyakov loop ([2]) in that it contains the scalar fields 
of the A/" = 4 theory as well as the gauge field. In 
that case, the boundary of the fundamental string world- 
sheet is located on the loop contour placed at the asymp- 
totic boundary of AdS^. In the zero temperature Yang- 
Mills theory defined on a spatial R'^, an interesting phe- 
nomenon occurs for loops in representations where the 
number of boxes k in the Young tableau is large so that 
jj: is finite in the large A^ limit. The dual fundamental 
string worldsheet is rep laced by a D-brane with world- 
volume electric fiux [l0l|-[T3|. This was found by studying 
highly supersymmetric i-BPS loops, where some results 
are known for all values of the coupling constant [iJl ■ For 
the anti-symmetric representation, the dual is a D5-brane 
whose world volume is a direct product of ^4^5*2 C AdS^ 
and C S^. For a symmetric representation, it is a 
D3-brane with world volume AdS2 x C AdS5. It is 
interesting to ask whether these D-branes exist in the 
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finite temperature geometry where they would be dual 
to a gauge theory loop linking periodic Euclidean time. 
This question has already been studied by Hartnoll and 
Kumar who looked for solutions of the appropriate 
Born-Infcld actions on the black hole background. For 
the D5-brane wrapped on S'^ C which corresponds 
to a totally antisymmetric representation on the gauge 
theory side, there seem to be solutions for any with 
the usual cutoff at k = N dictated by the maximum size 
of an antisymmetric representation on the gauge theory 
side and a maximum radius for embedding in on 
the supergravity side. However, in the case of the D3- 
brane which should correspond to a totally symmetric 
representation, Hartnoll and Kumar could not find any 
solutions at all. This fundamental difference between the 
two cases is what motivated our work on the gauge the- 
ory which we shall now summarize. Afterward, we will 
revisit the question of the supergravity side. We em- 
phasize that in supergravity we are studying the dual 
of the Maldacena- Wilson loop whereas in gauge theory 
our analysis is limited to the Polyakov loop. Both are 
governed by the center symmetry and both can become 
non-zero at the decnfinement transition. At high temper- 
ature, due to decoupling of the scalar fields, they should 
become similar. 

To study large totally symmetric or totally antisym- 
metric representations, it is convenient to obtain the 
characters in ([2|) from generating functions. 



^ (Tr5,/^, C/) = ^ dt- 



„TTrln(lTi(7)^ 



, (3) 



where the upper/lower sign is for symmet- 
ric/antisymmetric representations (Sk/Ak) respectively 
and the contour in the integral over t encircles the 
origin. In the large N limit, these integrals can be 
computed using two saddle point approximations. The 
first occurs while integrating over unitary matrices in 
([2]) . Because of the gauge symmetry, this is an eigenvalue 
model - the gauge symmetry can be used to diagonalize 



the eigenvalue density. 



^TJV/",d</'P(</')ln(lTte'*) 

27riife+i 



(4) 



The second use of a saddle-point approximation is to eval- 
uate the integral over t in (|4]). Let t satisfy the saddle- 
point equation 



1 T te*' 



k 
N 



(5) 



J —7T 

Then, the free energy is given by 

Ts./A, = ± dcj^pW ln(l T te'^) + ^ In £ . (6) 

The functions in ([5]) are related to the resol- 

vent of the matrix model and are holomorphic functions 
of t with cut singularities on the unit circle determined 
by the support of p{(j))- 

Let us begin with the symmetric representation Sk ■ We 
shall consider three examples of eigenvalue distributions. 
First, the confining phase has Pconf = ^^ki^) van- 
ishes if \t\ < 1 and is —1 if \t\ > 1. This is the expected 
discontinuity at the unit circle. Eq. ([5]) has solutions only 
when = 0, consistent with confinement. 

As a second example consider p{(j)) ~ 15^ (1 + 2pcos0). 
p = jf{Tr U) = J d4)p{4))e^'^ is the fundamental represen- 
tation loop. Positivity of the density requires < p < \. 
This distribution depends on cj) and therefore is decon- 
fined. While it is not realistic for = 4 Yang-Mills the- 
ory, it does occur in the strong-coupling phase of large 
N 2-dimensional lattice Yang-Mills theory p^ . 

There is one solution of Rs^{t) = in the region \t\ < 
I at i = /p. (If frackN and p are such that |t| > 1, 
both and F^^. should be extended there by analytic 
continuation.) The free energy is 



I = — in 



k/N 



ep 



(7) 



U 



diag[e* 



.,e' 



At large N, the eigenvalues where e = 2.718 . . .. Ts^ has the interesting feature that, 



become classical variables and their distribution is found 
by minimizing Scs phis a Jacobian from the unitary in- 
tegral measure. As long as fc << iV^, the loop operators 
in ([3]) do not modify the eigenvalue distribution which is 
given by a density pitp). p{(t))d(l) is times the number 
of eigenvalues between (f> and cj) -\- dcj) and is normalized, 
J_^d(l)p{(l)) = 1. Center symmetry is now an invariance 
under a simultaneous translation of all eigenvalues, 
(f'a — > (/)a+constant. In the center-symmetric confined 
phase, the distribution is translation invariant, eigen- 
values arc uniformly distributed on the unit circle and 
Pconf = In the de-confined phase, the eigenvalues are 
clumped. We will assume their distribution is symmetric 
about zero so that p{(j)) ~ p{—(j)). In the large TV limit 
the expectation values in Eq. ([3]) are computed using 



as is increased, it changes sign from negative to pos- 
itive. This results in a phase transition, which occurs 

5fc is 



N 



when 

negative and the loop expectation value. 



ep. When ^ < (^)^^^^, r, 



IS expo- 



nentially large. When > ^Sk is positive and 

the loop vanishes for iV — > 00. This phase transition im- 
plies that, even in the deconfined phase, sufficiently large 
symmetric representations are still confined. 

At this point, the reader might wonder how the expec- 
tation value of a unitary matrix can grow exponentially. 
The exponential comes from the traces needed to get the 
large representation and which give k and N dependent 
factors. Note that we did not normalize the loop (which 
would divide by a fc and TV-dependent factor). In Rcf.- 
p^ . it was shown that it is the un- normalized loop that 
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should be compared with supcrgravity. which is our even- 
tual aim. 

As a check of the saddle-point approximation to the 
t-integral in this simple example, observe that, if for the 
moment we assume that k and N are finite, we can in- 
tegrate (HI exphcitly to get e~^^^fc = ■^p'^- Using the 
Stirling formula and taking k N —f oo reproduces ([7]). 

We note that the presence of the phase transition is a 
universal property of the confining phase. From ^ and 
we see that ^^^J^^^ = Inf where i solves ([S]). Further, 

by inspecting ([6]) we see that Ts^. is real and negative 
when t is real and t < 1. As t increases, F^^. decreases to 
a minimum at i = 1, then begins increasing in the region 
t > 1 and eventually becomes positive, t increases with 
throughout this region. 

To see this behavior in another example, consider the 
semi-circle distribution which, for < 2 arcsin-y2 — 2p, 
is 



, ,/2 - 2w- sin^ 
7r(2 - 2p) V ^ 2 



(8) 



and which vanishes in the gap 2 arcsin ^2 — 2p < |0| < tt. 
We still use the fundamental loop, p, as a parameter and 
now i < J5 < 1. This is the distribution in the weak 
coupling phase of 2-dimensional lattice Yang-Mills the- 
ory p^ . It is also an approximation to the deconfined 
distribution for weakly coupled Af ~ 4 Yang-Mills theory 
For sufficiently weak coupling, it could be accu- 
rate near the phase transition where p = ^. The saddle 
point computation can be done explicitly near t=0 and 
analytically continued. The free energy is 



Ts^ = (26'cosh6'- sinh( 



sinh 6 + Vsinh^ 9 + 2 - 2p 



2-2p 



- In 



sinh 9 + v^sinh^ 9 + 2 - 2p 



2-2p 



where 9 is defined by t 
saddle-point equation 



,261 



and is determined by the 



1 



N 



cosh0 



sinh 9 + \/sinh^ 9 + 2-2j^ 
2-2p 



(10) 



which can be solved for sinh(0). The free energy is zero 
when fc = 0, negative for small fc, goes to zero at a critical 
and is positive thereafter. This is so for any value 
of p in the allowed range. A graph of F^^. versus for 
p = 0.51 is plotted in Fig. 1. With this value of p, the free 
energy becomes positive at ~ 0.50 which corresponds 
1.3. 

Now, we consider the antisymmetric representation. 
For a large class of distributions gapped around (f> = -k 

and with '^^'^'i > 0, which includes the semi-circle dis- 
tribution ([S]), we can argue that F^^ is always negative 
and the phase transition that we are discussing does not 



*° A'crit 



N 



FIG. 1: The free energy F^^ (S) as a function of in the 
semi-circle distribution with p = 0.51. 



occur. To begin, by changing variables in ([6]), we observe 
that F^j. = Tj\^_f,. This symmetry is reflected in the 
saddle-point equation ^ which, using our assumption 
that = p{—(f>), can be re- written as 



i- 7 — 

d(t)p{4>) — 



t 2 e 



f2e'2-|-i 26 



k 
N 



(11) 



and imphes i{k/N) = - k/N). The free energy 

d4>p{(j)) In (i^e 



+ 



hit 
(12) 

is symmetric under — > 1^ Moreover, with a gapped 
distribution, p(7r) ~ 0, and the integral in (fTTjl is contin- 



uous at i = 1. From 



At) 



> 0, i(fc) is monotone, and 
one can see in (fTT|) that £ = corresponds to = 0, 



since 



oo to 

dr 



N 

k 



= 1 and t = 1 to 



N 



N 

Furthermore, 



d^r. 



d{k/N) ~ lll^(^); d{k/Ny^ 



> 0, thus F_4j. is a convex 



function which decreases from to a negative minimum 
as goes from to ^ and then increases back to zero 
(9) when goes from i to 1. F^,, does not become positive 



and there is no phase transition of the kind that we found 
for symmetric representations. When the distribution is 

dRA^ (t) 



ungapped, or when 



dt 



becomes negative (for exam- 



ple when p < 0), interesting behavior can occur. We put 
off a discussion of it to a future investigation. 

We have found a difference between the symmetric 
and antisymmetric representation Polyakov loops which 
is qualitatively similar to the one found by HartnoU and 
Kumar [l5| for the dual objects in supcrgravity: the an- 
tisymmetric loop is non-zero in the deconfined phase for 
all allowed and the dual D5-brane exists whereas the 
gauge theory symmetric representation loop has a phase 
transition. The numerical search for the dual D3-brane in 
[Tst combined with analytic arguments at large k = ^ 
found no solution. If we take this to mean that the ex- 
pectation value of the gauge theory quantity vanishes, 
it suggests that, at strong coupling, the critical value 
goes to zero faster than 

We can also examine the alternative that the phase 
transion occurs for a value of k, so small that solutions 
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are missed in the numerical analysis. The gravity back- 
ground is the asymptotically AdS^ x black hole metric 



fir)df + 



+ {dx' + sin^ xd^ 



dni 



(13) 

where /(r) = {r^ — r'^){r'^ + 1 + r\)/r'^. The horizon 
is located at r = r-|_ and time is identified periodically 



with period 



2H +1 



the inverse Hawking temperature of 



the black hole. The D3-brane world volume embedded 
in this geometry is 



f{r)dt' + ^ + r^x'irfdr' + sin^ xir)dnl) 



(14) 

where t, r, [9, (f>) G are the world-volume coordinates 
which coincide with spacetime coordinates. The Born- 
Infeld action is S" = — f dtdrL, where 



L 



2 • 2 

r sm X 



47r2 

l + r\f{r){xr-—F''-r\' 



(15) 

The first term is y^det{g + F), with gab the worldvolume 
metric and Fab is the worldvolume electromagnetic field. 
The last term is the Chern-Simons term. Consistent with 
symmetries and equations of motion, one can take Fat 
with one nonzero component Frt = A'^ (r) . The canonical 
momentum 11 = dL/dA^{r) is a constant equal to the 
number of units of electric flux, k. Solving for F and 
substituting into the equation of motion for x yields 



„ 4 . 3 ^l + rVi^f 

2r sm xcosx- 



•\/7^iI?X■ 



4r sin X 




r^Iix'f 



(16) 



We fix the boundary condition at r — > oo to match the 
zero temperature i-BPS D3-brane. It is a solution of 
the same equation with fir) = and where, to get 
the Poincare coordinates, (sin x, cos x) are replaced by 
(x, !)■ Then, xi''') = 7 is an exact solution of ((TBI) ^^^^ 
the brane geometry is a simple direct product of AdS2 
with radius vT+^ and S*^ with radius k [i3| . We seek 
solutions with the asymptotic behavior xi^") ~ 7 for large 
r. By studying the large r regime, it is easy to see that 
there is no solution of (jl6p for x which goes to zero at 
least as fast as r"^ unless n is non-zero. By studying 
the region near the horizon, we can see that there is no 
solution in the large k limit. So, if there is a solution at 
all, it will only exist if kappa is non-zero but not large. 
We have attempted to solve (fT6|) numericaly with small 
values of k. We have positive evidence for a solution in 
a corner of the parameter space obtained by taking the 
infinite temperature limit (replacing r with rL, r+ by 
r+L, X with x/L and taking L 00). The resulting 
differential equation has an exact solution for ?■+ ~ 0, 
X = K/{r + b) + 0{r'^) where b is an integration constant. 
Restoring r+ > 0, we employed a shooting technique to 
look for solutions which asymptote to x- With k = 0.001, 
there appears to be a solution at ~ 10^°. Our work is 
on-going and we shall present the details elsewhere. 
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